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INTRODUCTION 

Aquatic  invertebrates  are  aptly  applied  to  bioassessment  since  they  are  known 
to  be  important  indicators  of  stream  ecosystem  health  (H)Ties  1970).  Long  lives, 
complex  life  cycles  and  limited  mobility  mean  that  there  is  ample  time  for  the  benthic 
community  to  respond  to  cumulative  effects  of  environmental  perturbations.  This  report 
summarizes  data  collected  in  August  -  September  2003  from  9  sites  in  the  Madison 
River  drainage.   Sites  were  located  in  Madison  and  Beaverhead  Counties,  Montana  and 
lie  within  the  Montana  Valley  and  FoothiU  Prairies  and  Northern  Rockies  ecoregions 
(Woods  et  al.  1999). 

A  mvdtimetric  approach  to  bioassessment  such  as  the  one  applied  in  this  study 
uses  attributes  of  the  assemblage  in  an  integrated  way  to  measure  biotic  health.  A 
stream  with  good  biotic  health  is  "...a  balanced,  integrated,  adaptive  system  having  the 
full  range  of  elements  and  processes  that  are  expected  in  the  region's  natural 
environment..."  (Karr  and  Chu  1999).  The  approach  designed  by  Plafkin  et  al.  (1989) 
and  ad^ted  for  use  in  the  State  of  Montana  has  been  defined  as  ". . .  an  array  of 
measures  or  metrics  that  individually  provide  information  on  diverse  biological 
attributes,  and  when  integrated,  provide  an  overall  indication  of  biological  condition." 
(Barbour  et  al.  1995).  Community  attributes  that  can  contribute  meaningfully  to 
interpretation  of  benthic  data  include  assemblage  structure,  sensitivity  of  community 
members  to  stress  or  pollution,  and  functional  traits.  Each  metric  component 
contributes  an  independent  measure  of  the  biotic  integrity  of  a  stream  site;  combining 
the  components  into  a  total  score  reduces  variance  and  increases  precision  of  the 
assessment  (Fore  et  al.  1996).  Effectiveness  of  the  integrated  metrics  depends  on  the 
applicability  of  the  underlying  model,  which  rests  on  a  fovindation  of  three  essential 
elements  (BoUman  1998a).  The  first  of  these  is  an  appropriate  stratification  or 
classification  of  stream  sites,  typically  by  ecoregion.  Second,  metrics  must  be  selected 
based  upon  their  ability  to  accurately  express  biological  condition.  Third,  an  adequate 
assessment  of  habitat  conditions  at  each  site  to  be  studied  enhances  the  interpretation  i 

of  metric  outcomes. 

Implicit  in  the  multimetric  method  and  its  associated  habitat  assessment  is  an  \ 

assumption  of  correlative  relationships  between  habitat  measures  and  the  biotic  ,^ 

metrics,  in  the  absence  of  water  quality  impairment.  These  relationships  may  vary 
regionally,  requiring  an  examination  of  habitat  assessm.ent  elements  and  biotic  metrics 
and  a  test  of  the  presumed  relationship  between  them.  BoUman  ( 1998a)  has  studied  the 
assemblages  of  the  Montana  Valley  and  Foothill  ftairies  (MVFF^  ecoregion  and  has  j 

reconmiended  a  battery  of  metrics  appficable  to  the  montane  ecoregions  of  western  j- 

Montana.  This  metric  battery  has  been  shown  to  be  sensitive  to  impairment,  related  to 
measures  of  habitat  int^rity,  and  consistent  over  replicated  samples.  However,  scoring 
criteria  developed  for  the  MVFP  ecoregion  metric  battery  may  not  be  appropriately 
sensitive  for  streams  of  the  Northern  Rockies  ecoregion.  Additional  work  may  be  needed 
to  calibrate  scoring  to  montane  regions,  thus  impairment  classifications  must  be 
interpreted  with  care.  For  comparison,  the  bioassessment  index  foiond  in  the  Montana 
Department  of  Environmental  Quality  (Montana  DEQ)  Standard  Operating  Procedures  /  i 

(Bukantis  1998)  is  also  scored  and  presented.  However,  this  index  has  not  been  tested 
for  sensitivity  to  impairment,  relationship  to  measures  of  habitat  integrity,  or 
consistency  over  replicated  samples. 

In  this  report,  3  assessment  methods  are  used:  first,  taxonomic  data  is 
translated  in  to  2  bioassessment  indices,  and  metric  scores  in  each  index  are  summed 
to  derive  impairment  classifications  and  use  support  designations.  Second,  a  narrative 
interpretation,  based  on  the  author's  professional  judgment  is  given.  Metric 
performance  and  taxonomic  data  are  both  applied  to  this  analysis.  Third,  the  model  of 
Barbour  and  Stribling  (1991)  is  applied  to  bioassessment  and  habitat  assessment 

scores.  While  the  interdependence  of  these  methods  is  obvious,  since  the  same  data  are      

used  for  all,  some  degree  of  independence  is  maintained  throughout  the  analysis. 
Narrative  interpretations  are  given  without  regard  to  the  bioassessment  index  result 


and  without  reference  to  habitat  assessment.  Since  indices  are  summations,  they  can 
often  mask  evidence  of  impairment;  the  narratives  attempt  to  expose  the  p>otential 
shortcomings  of  the  indices.  Similarly,  graphing  the  association  between  habitat 
assessment  scores  and  bioassessment  scores  using  the  model  of  Barbour  and  Stribling 
can  provide  clues  to  offer  support  or  possible  refutation  of  the  conclusions  of  the 
narrative  analysis. 

METHODS 

Samples  were  collected  in  August  -  September  2003  by  personnel  of  the 
Montana  Department  of  Environmental  Quality  (Montana  DEQ).  Sample  designations 
and  site  locations  are  indicated  in  Table  1  and  on  Figure  1 .  The  site  selection  and 
kicknet  sampling  method  employed  were  those  recommended  in  the  Montana  DEQ 
Standard  Operating  Procedures  for  Aquatic  Macroinvertebrate  Sampling  (Bukantis 
1998).  Aquatic  invertebrate  samples  were  delivered  to  Rhithron  Associates,  Inc., 
Missoula,  Montana,  for  laboratory  and  data  analyses.  In  the  laboratory,  the  Montana 
DEQ-recommended  sorting  method  was  used  to  obtain  subsamples  of  at  least  300 
organisms  from  each  sample,  \dien  possible.  Organisms  were  identified  to  the  lowest 
possible  taxonomic  levels  consistent  with  Montana  DEQ  protocols. 

Table  1.  Sample  designations  and  locations.  Sites  are  listed  in  upstream-to- 
downstream  order.  Madison  River  drainage.  August  -  September  2003. 


Site  ID 


Statloo  ID  Activity  ID  Location  Description 


Collection 
Date 


Latitude 


Longltade 


BUFUC 
01 

ELKR 
01 

ELKR 
02 

autlc 

01 

OAZLC 
01 

OAZLC 
02 

BLIfSC 
03 

Bunc 

Ol 

BUfSC 
04 


M06BUFDCX>I       03-U345-M 


M06ELKR01 
M06ELKR02 
M06ANTLC01 
M06QAZLC01 


03-U341-M 
03-U342-M 
03-U340-M 
03-U343-M 


M06GAZLC02   03-U344-M 


M06BLNSC03   03-U355-M 


M06BLNSC01   03-U353-M 


M06BLNSC04        03-U356-M 


BUFORD  CREEK  1/8  MI 
U/S  FROM  MOUTH  OF 
W  FK  MADISON 

ELK  RTVER  4  MI  UP 
FSRD  1209 

ELK  RIVER  @  END  OF  W 
FK  MADISON  R  RD 

ANTELOPE  CREEK 
UPPER  REACH  ON 
ANTELOPE  BASIN  RD 
GAZELLE  CREEK  1/3  MI 
D/S  OF  XING  GAZELLE 
CR  AND  FSRD  1200 
GAZELLE  CREEK  200 
YDS  U/S  FROM  FS  209 
XlNQ  @  GAZELLE  CR 
BLAINE  SPRING  CREEK 
BELOW  FISH 
HATCHERY 

BLAINE  SPRING  CREEK 
1/3  MI  U/S  FROM 
MOUTH 

BLAINE  SPRING  CREEK 
@  VARNEY 


8/28/03 

44-45- 18 

111-51-27 

8/26/03 

44-48-0 

111-39-51 

8/26/03 

44-47-28 

111-30-4 

8/25/03 

44-41-3 

111-31-40 

8/27/03 

44-52-8 

111-41-15 

8/27/03 

44-53-8 

111-35-18 

9/11/03 

45-12-55 

111-47-30 

9/10/03 

45-17-6 

111-45-26 

9/11/03 

45-18-59 

111-45-25 

Figure  1.  Approximate  sampling  locations.  Madison  River  drainage.  August 

September  2003 
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Two  bioassessment  indices  were  employed  in  evaluating  the  data.  The  first  index 
employed  was  developed  for  streams  of  western  Montana  ecoregions  (Bollman  1998a). 
Each  metric  used  in  this  index  was  tested  for  its  response  or  sensitivity  to  varying 
degrees  of  hiunan  influence.  Correlations  have  been  demonstrated  between  the  metrics 
and  various  symptoms  of  human-caused  impairment  as  expressed  in  water  quality 
parameters  or  instream,  streambank,  and  stream  reach  morpholc^c  features.  Metrics 
were  screened  to  minimize  variability  over  natural  environmental  gradients,  such  as  site 
elevation  or  sampling  season,  viaich  might  confound  interpretation  of  results  (Bollman 
1998a).  The  multimetric  index  incorporates  multiple  attributes  of  the  sampled 
assemblage  into  an  int^rated  score  that  accurately  describes  the  benthic  community  of 
each  site  in  terms  of  its  biologic  int^rity.  In  addition  to  the  metrics  included  in  the 
index,  other  metrics  shown  to  be  applicable  to  biomonitoring  in  other  regions  (Kleindl 
1995,  Patterson  1996,  Rossano  1995)  were  used  for  descriptive  interpretation  of  results. 
These  metrics  include  the  number  of  "dinger"  taxa,  long-lived  taxa  richness,  the  percent 
of  predatory  organism.s,  and  others.  They  are  not  included  in  the  int^rated 
bioassessment  score,  however,  since  their  performance  ia  western  Montana  ecoregions 
is  unknown.  However,  the  relationship  of  these  metrics  to  habitat  conditions  is  intuitive 
and  reasonable. 

The  six  metrics  included  in  the  bioassessment  index  used  for  MVFP  sites  in  this 
study  were  selected  because,  both  individually  and  as  an  integrated  metric  batteiy,  they 
are  robust  at  distingviishing  impaired  sites  from  relatively  unimpaired  sites  (Bollman 
1998a).  In  addition,  they  are  relevant  to  the  kinds  of  impacts  that  are  present  in  the 
Madison  River  drainage.  They  have  been  demonstrated  to  be  more  variable  with 
anthropogenic  disturbance  than  with  natural  environmental  gradients  (Bollman  1998a). 
Each  of  the  six  metrics  developed  and  tested  for  western  Montana  ecoregions  is 
described  below. 

1.  Bphemeroptera  (mayfly)  taxa  richness.  The  number  of  mayfly  taxa 
declines  as  water  quality  diminishes.  Impairments  to  water  quality  which  have  been 
demionstrated  to  adversely  affect  the  ability  of  mayflies  to  flourish  include  elevated 
water  tetajjeratures,  heavy  metal  contamination,  increased  turbidity,  low  or  high 
pH,  elevated  specific  conductance  and  toxic  chemicals.  Few  m.ayffy  species  are  able 
to  tolerate  certain  disturbances  to  instream  habitat,  such  as  excessive  sediment 
deposition. 

2.  Plecoptera  (stonefly)  taxa  richness.  Stoneflies  are  particvdarly 
susceptible  to  impairments  that  afiiect  a  stream  on  a  reach-level  scale,  such  as  loss 
of  riparian  canopy,  streambank  instability,  channelization,  and  alteration  of 
morphological  features  such  as  pool  firequency  and  function,  riffle  development  and 
sinuosity.  Just  as  all  benthic  organisms,  they  are  also  susceptible  to  smaller  scale 
habitat  loss,  such  as  by  sediment  deposition,  loss  of  interstitial  spaces  betvpeen 
substrate  particles,  or  unstable  substrate. 

3.  Trlchoptera  (caddisily)  taxa  richness.  Caddisfly  taxa  richness  has  been 
shown  to  decline  when  sediment  deposition  affects  habitat.  In  addition,  the  presence 
of  certain  case-building  caddisflies  can  indicate  good  retention  of  woody  debris  and 
lack  of  scouring  flow  conditions. 

4.  Number  of  sensitive  taxa.  Sensitive  taxa  are  generally  the  first  to 
disappear  as  anthropogenic  disturbances  increase.  The  list  of  sensitive  taxa  used 
here  includes  oiganisms  sensitive  to  a  wide  range  of  disturbances,  including  warmer 
water  temperatures,  organic  or  nutrient  pollution,  toxic  pollution,  sediment 
deposition,  substrate  instabilityr  and  others.  Unimpaired  streams  of  western 
Montana  typically  support  at  least  four  sensitive  taxa  (Bollman  1998^. 

5.  Percent  filter  feeders.   FUter-feediog  organisms  are  a  diverse  group; 
-    they  capture  small  particles  of  organic  matter,  or  organically  enriched  sediment 

material,  fi'om  the  water  column  by  means  of  a  variety  of  adaptations,  such  as 


silken  nets  or  hairy  appendages.  In  forested  montane  streams,  filterers  are  expected 
to  occur  in  insignificant  nxunbers.  Their  abvmdance  increases  when  canopy  cover  is 
lost  and  \rtien  water  temperatxires  increase  and  the  accompanying  growth  of 
filamentous  algae  occurs.  Some  filtering  organisms,  specifically  the  Arctopsychid 
caddisfiies  (Arctopsyche  spp.  and  Parapsyche  spp.)  build  silken  nets  with  large  mesh 
sizes  that  capture  small  organisms  such  as  chironomids  and  early-instar  mayflies. 
Here  they  are  considered  predators,  and,  in  this  study,  their  abundance  does  not 
contribute  to  the  percent  filter  feeders  metric. 

6.    Percent  tolerant  taxa.  Tolerant  taxa  are  ubiquitous  in  stream  sites,  but 
when  disturbance  increases,  their  abundance  increases  proportionately.  The  list  of 
taxa  used  here  includes  organisms  tolerant  of  a  wide  range  of  disturbances, 
including  warmer  water  temperatures,  organic  or  nutrient  pollution,  toxic  pollution, 
sediment  deposition,  substrate  instability  and  others. 

A  second  index  was  appKed;  depending  on  the  site  location,  the  Montana  DEQ 
index  for  mountain  streams  or  for  foothill  streams  was  used.  These  are  the  metric 
batteries  recommended  in  the  Montana  DEQ  standard  operating  procedures  (Bukantis 
1998).  Both  of  these  indices  should  be  considered  provisional,  since  correlative 
relationships  between  them  and  meaningfiil  measures  of  habitat  condition  and  water 
quality  have  not  been  evaluated.  Thus,  conclusions  concerning  bioassessment  based 
upon  these  metrics  must  be  regarded  as  tentative. 

Scoring  criteria  for  each  of  the  metrics  in  the  indices  applied  in  this  study  are 
presented  in  Tables  2a,  2b,  and  2c.  Metrics  differ  in  their  possible  value  ranges  as  well 
as  in  the  direction  the  values  move  as  biological  conditions  change.  For  example, 
Ephemeroptera  richness  values  may  range  fi-om  zero  to  ten  taxa  or  higher.  Larger 
values  generally  indicate  favorable  biotic  conditions.  On  the  other  hand,  the  percent 
filterers  metric  may  range  firom  0%  to  100%;  in  this  case,  larger  values  are  negative 
indicators  of  biotic  health.  To  facilitate  scoring,  therefore,  metric  values  were 
transformed  into  a  single  scale.  The  range  of  each  metric  has  been  divided  into  four 
parts  and  assigned  a  point  score  between  zero  and  three.  A  score  of  three  indicates  a 
metric  value  similar  to  one  characteristic  of  a  non-impaired  condition.  A  score  of  zero 
indicates  strong  deviation  fi-om  non-impaired  condition  and  suggests  severe  degradation 
of  biotic  health.  Scores  for  each  metric  were  summed  to  give  an  overall  score,  the  total 
bioassessment  score,  for  each  site  in  each  sampling  event.  These  scores  were  expressed 
as  the  percent  of  the  maximum  possible  score.  The  total  bioassessment  score  for  each 
site  was  expressed  in  terms  of  use-support.  Criteria  for  use-support  designations  were 
developed  by  Montana  DEQ  and  are  presented  in  Table  3a.  Scores  were  also  translated 
into  impairment  classifications  according  to  criteria  outlined  in  Table  3b. 

In  this  report,  certain  other  metrics  were  used  as  descriptors  of  the  benthic 
community  response  to  habitat  or  water  quality  but  were  not  incorporated  into  the 
bioassessment  metric  batteiy,  either  because  they  have  not  yet  been  tested  for 
reliability  in  streams  of  western  Montana,  or  because  results  of  such  testing  did  not 
show  them  to  be  robust  at  distinguishing  impairment,  or  because  they  did  not  meet 
other  requirements  for  inclusion  in  the  metric  battery.  These  metrics  and  their  use  in 
predicting  the  causes  of  impairment  or  in  describing  its  effects  on  the  biotic  community 
are  described  below. 


Table  2a.  Metrics  and  scoiing  criteria  for  bioassessment  of  montane  streams  of  Western 
Montana  (Bollman  1998£^. 

Score 
Metric  3  2  1  0 

Ephemeroptera  taxa  richness 
Plecopterm  taza  richness 
Trichopterm  taxa  richness 
Sensitive  tmxa.  richness 
Percent  fllterers 
-  ■'"■    Percent  tolerant  taxa 


>5 

5-4 

3-2 

<2 

>3 

3-2 

1 

0 

>4 

4-3 

a 

<2 

>3 

3-2 

1 

0 

0-5 

5.01  -  10 

10.01-25 

>25 

0-5 

5.01  -  10 

10.01-35 

>35 

Table  2b.  Provisional  metrics  and  scoring  criteria  for  bioassessment  of  streams  of  the 
Mountain  Ecoregions  (Bukantis  1998). 


Score 

Metric 

3 

a 

t 

0 

Taza  richness 

>28 

28 

-24 

24 

-  19 

<  19 

EPT  richness 

>  19 

19 

-  17 

17 

-  15 

<  15 

Blotic  index 

<3 

3 

-4 

4- 

-5 

>S 

%  dominant  taxon 

<25 

25 

-  35 

35 

-45 

>45 

%  collector 

<60 

60 

■70 

70 

-80 

>80 

%  scraper  + 
shredder 

>55 

55- 

■40 

40 

-25 

<25 

%EPT 

>70 

70- 

■55 

55 

-40 

<40 

Table  2c.  Provisional  metrics  and  scoring  criteria  for  bioassessment  of  streams  of  the 
Montana  Valley  and  Foothill  Prairies  Ecoregion  (Bukantis  1998). 


Score 

Metric 

3 

a 

1 

0 

Taxa  richness 

>28 

28 

-  21 

20 

■  14 

<  14 

EPT  richness 

>  14 

14 

-  13 

12 

■  11 

<  11 

Biotlc  index 

•<4 

4 

-5 

5- 

-6 

>6 

%  dominant  taxon 

<30 

30 

■40 

40 

•50 

>50 

%  collector 

<60 

60 

■75 

75 

■90 

>90 

%  scraper  * 
shredder 

>30 

30 

-20 

20 

■  10 

<  10 

%EPT 

>60 

60 

-45 

45- 

■  30 

<30 

Percent 

Hydropsychidae 

of 

<75 

75- 

-85 

85  ■ 

95 

>95 

Trichoptera 

TableSa.  Criteria  for  the  assignment  of  use-support  classifications  /  standards  violation 
thresholds  (Bukantis  1998). 

%  Comparability  to  reference 

Use  support 

>75 

25-75 

<25 

Full  support— standards  not  violated 

Partial  support— moderate  impairment— standards  violated 

Non-support— severe  impairment— standards  violated 

TableSb.  Criteria  for  the  assigimient  of  impairment  classifications  (Plafkin  et  al.  1989). 

%  Comparability  to  reference 

Classification 

>83 

54-79 
21-50 
<17 

nonimpaired 
slightly  impaired 
moderately  impaired 
severely  impaired 

The  modified  biotic  index.  This  metric  is  an  adaptation  of  the  Hilsenhofi"  Biotic 
Index  (HBI,  HUsenhoff  1987),  which  was  originally  designed  to  indicate  oiganic 
enrichment  of  waters.  Values  of  this  metric  are  lowest  in  least  impacted 
conditions.  Tzixa  tolerant  to  saprobic  conditions  are  also  generally  tolerant  of 
warm  water,  fine  sediment  and  heavy  filamentous  algae  growth  (Bollman 
1998b).  Loss  of  canopy  cover  is  often  a  contributor  to  higher  biotic  index  values. 
The  taxa  values  used  in  this  report  are  modified  to  reflect  habitat  and  water 
quality  conditions  in  Montana  (Biikantis  1998).  Ordination  studies  of  the 
benthic  faiuia  of  Montana's  foothill  prairie  streams  showed  that  there  is  a 
correlation  between  modified  biotic  index  values  and  water  temperature, 
substrate  embeddedness,  and  fine  sediment  (Bollman  1998a).  In  a  study  of 
reference  streams,  the  average  value  of  the  modified  biotic  index  in  least- 
impaired  streams  of  western  Montana  was  2.5  (Wisseman  1992). 
Taxa  richness.  This  metric  is  a  simple  count  of  the  number  of  unique  taxa 
present  in  a  saiaple.  Average  taxa  richness  in  samples  from  reference  streams  in 
western  Montana  was  28  (Wisseman  1992).  Taxa  richness  is  an  expression  of 
biodiversity,  and  generally  decreases  with  degraded  habitat  or  diminished  water 
quality.  However,  taxa  richness  may  show  a  paradoxical  increase  vt^en  mild 
nutrient  enrichment  occurs  in  previously  oLlgotrophic  waters,  so  this  metric 
must  be  interpreted  with  caution. 

Percent  predators.  Aquatic  invertebrate  predators  depend  on  a  reliable  source  of 
invertebrate  prey,  and  their  abundance  provides  a  measure  of  the  trophic 
complexity  supported  by  a  site.  Less  disturbed  sites  have  more  plentiful  habitat 
niches  to  support  diverse  prey  species,  which  in  turn  support  abundant 
predator  species. 

Number  of  "dinger"  taxa.  So-called  "dinger"  taxa  have  physical  adaptations  that 
allow  them  to  cling  to  smooth  substrates  in  rsqjidly  flowing  water.  Aquatic 
invertebrate  "dingers"  are  sensitive  to  fine  sediments  that  fill  interstices  between 
substrate  particles  and  eliminate  habitat  complexity.  Animals  that  occupy  the 
hyporheic  zones  are  included  in  this  group  of  taxa.  Expected  "dinger^  taxa 
richness  in  unimpaired  streams  of  western  Montana  is  at  least  14  (Bollman 
1998b). 

Ntunber  of  long-lived  taxa.  Long-lived  or  semivoltine  taxa  require  more  than  a 
year  to  completely  develop,  and  their  numbers  decline  wiien  habitat  emd/or 
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water  quality  conditions  are  unstable.  They  may  completely  disappear  if 
channels  are  dewatered  or  if  there  are  periodic  water  temperature  elevations  or 
other  interruptions  to  their  Ufe  cycles.  Western  Montana  streams  with  stable 
habitat  conditions  are  e^qiected  to  support  six  or  more  long-lived  taxa  (BoUman 
1998b). 


RESULTS 


Bioassessment 


Figure  2  summarizes  bioassessment  scores  for  aquatic  invertebrate  communities 
sampled  from  sites  in  the  Madison  River  drainage.  Tables  4a  and  4b  itemize  each 
contributing  metric  and  show  individual  metric  scores  for  each  site.  Tables  3a  and  3b 
above  show  criteria  for  use-support  categories  recommended  by  Montana  DEQ 
(Bukantis  1998)  and  impairment  classifications  (Plafkin  et  al.  1989) 

When  the  revised  assessment  method  is  applied  to  these  data,  scores  indicate 
that  4  of  the  9  sites  studied  were  unimpaired  and  fully  supported  designated  uses. 
These  wrere  ELKROl,  ANTLCOl,  GAZLCOl,  and  GAZLC02.  One  site,  ELKR02,  was 
slightly  impaired  but  fully  supportive  of  uses.  Site  BURDCOl  was  slightly  impaired  and 
only  partly  supported  designated  uses.  Moderate  impairment  and  partial  use  support 
was  indicated  at  2  sites:  BLNSC03  and  BLNSC04.  Site  BLNSCOl  was  severely  impaired 
and  did  not  support  uses. 

The  2  bioassessment  methods  did  not  produce  comparable  results  in  several 
cases;  differences  were  remarkable  at  several  sites,  notably  ELKR02  and  BLNSCOl. 
Sites  were,  however,  ranked  similarly,  with  the  exception  of  ELKR02  and  BLNSCOl. 

Fi£are  2.  Comparison  of  total  bioassessment  scores  (reported  as  percent  of  maximum 
score)  for  9  sites  in  the  Madison  River  drainage.  August  -  September  2003.  The  revised 
bioassessment  method  (BoUman  1998a)  was  used  to  calculate  scores  represented  by 
orange  bars,  and  the  Montana  DEQ  bioassessment  method  (Bukantis  1998)  was  used 
to  calculate  scores  represented  by  purple  bars. 
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Aquatic  invertebrate  communities 

Interpretations  of  biotic  integrity  in  this  report  are  made  without  reference  to 
results  of  habitat  assessments,  or  any  other  information  about  the  sites  or  watersheds 
that  may  have  accompanied  the  invertebrate  samples.  Interpretations  are  based  entirely 
on:  the  taxonomic  and  functional  composition  of  the  sampled  invertebrate  assemblages; 
the  sensitivities,  tolerances,  physiology,  and  habitus  information  for  iadividual  taxa 
gleaned  from  the  vmter's  research;  the  published  literature,  and  other  expert  sources; 
and  on  the  performance  of  bioassessment  metrics,  described  earlier  in  the  report,  which 
have  been  demonstrated  to  be  useful  tools  for  interpreting  potential  implications  of 
benthic  invertebrate  assemblage  composition. 

At  the  sampled  site  on  Buford  Creek  (BURDCOl),  fewer  ma5^es  than  expected 
were  present  in  the  sample,  and  the  biotic  index  value  (4.08)  was  higher  than  expected. 
These  findings  suggest  that  water  quality  may  have  been  degraded  by  nutrient  v 

enrichment,  elevated  water  temperatures,  or  both.  A  single  cold-stenotherm  taxon,  the 
midge  Cricotopus  nostocodadius,  was  collected.  Thirty-one  percent  of  the  collected 
invertebrates  were  midges.  The  abvmdance  of  the  caddisfly  Hydroptila  sp.  su^ests  the 
presence  of  filamentous  algae,  often  an  indicator  of  nutrient  enrichment.  Thus,  it  is 
possible  to  discern  evidence  of  both  types  of  water  quality  impairment  in  the 
invertebrate  community  data.  The  psychodid  fly  Pericoma  sp.  was  common  at  this  site, 
suggesting  carbonate  topology  in  the  watershed. 

Although  18  "dinger"  taxa  were  taken,  there  were  only  3  caddisfly  taxa  in  the 
sampled  assemblage.  It  is  not  possible  to  rule  out  disruption  of  instream  habitats  by 
sediment  deposition,  given  these  findings.  Since  14  predator  taxa  were  among  the  rich 
fauna  (43  taxa),  it  seems  likely  that  instream  habitats  of  diverse  types  were  available. 
Reach-scale  habitat  features,  such  as  channel  morphological  elements,  streambanks, 
and  riparian  zones,  may  have  been  disturbed  to  some  extent,  since  only  3  stonefly  taxa 
were  present.  The  site  supported  at  least  5  semivoltine  taxa;  surface  flow  must  have 
persisted  year-round,  and  recent  scouring  sediment  pulses  are  not  indicated.  All 
expected  functional  comiponents  were  present.  Notable  is  the  abundance  of  shredders, 
indicating  ample  riparian  inputs  of  large  organic  material,  and  hydrologic  conditions 
suitable  for  its  retention. 

The  upper  site  on  the  Elk  River  (ELKROl)  supported  a  diverse,  sensitive,  and 
functional  assemblage  characteristic  of  minimaUy  impaired  montane  sites.  The  low 
biotic  index  value  (1.93)  and  rich  mayfly  fauna  (12  taxa)  suggests  excellent  water 
quality.  The  presence  of  10  cold  stenotherm  taxsi,  including  Baetis  bicaudatus, 
Caudatella  sp.,  and  Deuterophlebia  sp.,  indicates  cold  temperatures. 

Six  caddisfly  taxa  and  21  "dinger"  taxa  were  coixnted  in  the  sample;  sediment 
deposition  did  not  likely  contaminate  stony  benthic  habitats  here.  Diverse  undisturbed 
instream  habitats  appear  to  be  indicated  by  the  rich  fauna;  thirty-eight  taxa  were        ? 
collected,  and  among  these,  12  were  predators.  Reach  scale  habitat  features  were  likely 
intact,  since  there  were  6  stonefly  taxa  in  the  assemblage.  Year-round  surface  flow  and 
the  lack  of  any  recent  catastrophes  were  indicated  by  the  5  long-lived  taxa  supported  at 
the  site.  Although  all  expected  functional  components  were  present,  the  abvindance  of 
scrapers  and  relative  paucity  of  shredders  may  imply  a  lack  of  riparian  shading  and 
allochthonous  inputs. 

The  downstream  site  on  the  Elk  River  (ELKR02)  supported  7  mayfly  taxa, 
suggesting  good  water  quality,  but  the  biotic  index  value  calcxdated  for  the  assemblage 
as  a  whole  was  4.28,  higher  than  expected,  and  su^estive  of  impaired  water  quality. 
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The  sampled  assemblage  was  dominated  (54%)  by  the  blackfly  SimuUum  sp.,  the 
abundance  of  which  has  apparently  skewed  the  biotic  index  value.  Blackfly  larvae  in 
this  genus  are  gregarious,  and  samples  can  be  swamped  by  them  when  a  cluster  is 
encovuitered  by  the  sampler.  When  the  calculation  is  made  without  including  Simuhum 
sp.,  the  resiilting  value  is  2.25.  Four  sensitive  cold-stenotherm  taxa,  including  Drunella 
doddsi  and  Caudatella  sp.  were  common  in  the  sample.  It  appears  that  water  quality 
was  probably  good,  and  water  temperatures  cold  at  this  site. 

Seventeen  "dinger'  taxa  and  8  caddisfly  taxa  were  taken  in  the  sample, 
suggesting  clean  stony  substrates  without  appreciable  sediment  deposition.  Other 
instream  habitats  were  likely  diverse  and  available,  since  at  least  29  distinct  taxa  were 
supported  at  the  site;  nine  of  these  taxa  were  predators.  The  stonefly  fauna  was  limited 
to  2  taxa.  This  may  be  associated  with  disturbance  to  riparian  vegetation,  streambank 
instability,  or  alterations  to  natural  chaxmel  morpholc^y.  Five  semivoltine  taxa  were 
counted,  implying  year-round  surface  flow.  All  expected  functional  components  vrere 
present.  The  dominance  of  SimuKum.  sp.  has  resvdted  in  the  overwhelming  abundance  of 
filter-feeders  among  the  feeding  groups  represented  here. 

The  site  on  Antelope  Creek  (ANTLCO 1)  supported  a  sensitive  and  diverse  benthic 
assemblage.  High  ma5rfly  taxa  richness  (7)  and  a  low  biotic  index  value  (2.82)  are 
suggestive  of  excellent  water  quality.  No  fewer  than  4  cold-stenotherm  taxa  were  present 
at  the  site,  including  the  dipteran  Glutops  sp.  and  the  caddisfly  Apatania  sp.  The 
turbeUarian  Polycelis  corvnata  was  common,  siiggesting  that  groixndwater  inputs 
augmented  flow  here. 

Six  caddisfly  taxa  and  14  "dinger'  taxa  vrere  collected.  It  is  unlikely  that 
sediment  deposition  interfered  with  the  accessibflity  of  stony  substrate  habitats.  The 
rich  fauna  (35  taxa)  and  numierous  predatory  taxa  (9)  imply  that  instream  habitats  of 
diverse  types  were  available  and  undisturbed.  There  were  fewer  stonefly  taxa  (3)  in  the 
sample  than  expected;  low  richness  in  this  order  may  be  associated  with  disturbance  to 
riparian  zones,  streambanks,  or  natxxral  channel  m.orphology.  Only  2  long-lived  taxa 
were  collected;  the  possibility  of  recent  dewatering,  scouring  sediment  pulse,  or  other 
catastrophic  event  cannot  be  ruled  out.  All  expected  functional  components  veere 
present;  the  abundance  of  collectors  is  related  to  the  dominance  of  the  elmid  beetle 
HeterHmmus  sp.  in  the  sample.  Since  this  animal  is  gregarious  in  habit,  their  numbers 
may  be  overrepresented  in  kick  net  collections. 

Both  sites  on  Gazelle  Creek  (GAZLCO 1  and  GAZLC02)  supported  sensitive, 
diverse,  and  functional  assemblages  characteristic  of  minimally  disturbed  montane 
streams.  Thirteen  taxa  collected  at  the  upper  site  (GAZLCO  1)  were  cold-stenotherms; 
these  included  the  nemourid  stonefly  Visoka  cataractae  and  the  caddisfly  Parapsyche 
elsis.  The  biotic  index  value  (2.46)  was  low,  and  the  mayfly  taxa  richness  (1 1)  high;  both 
of  these  findings  suggest  excellent  water  quality.  Polycelis  coronata  was  abundant, 
sviggesting  the  influence  of  groundwater  inputs  on  stream  flow. 

Twenty-four  "dinger"  taxa  and  1 1  caddisfly  taxa  vrere  among  the  collected 
assemblage,  making  it  seem  vety  unlikely  that  sediment  deposition  compromised  stony 
substrate  habitats.  In  all,  53  discreet  taxa  were  counted  in  the  sample,  and  22  of  these 
were  predatory.  Complex,  diverse  instream  habitats  are  suggested  by  the  richness  of  the 
fauna.  Reach-scale  habitat  features  were  likely  intact,  since  the  site  supported  at  least 
9  stonefly  taxa.  Dewatering  or  other  catastrophic  interruptions  to  long  life  cycles  were 
probably  not  among  recent  events  here;  four  semivoltine  taxa  were  present.  All  expected 
functional  components  were  coUected  in  appropriate  proportions. 

At  the  downstream  site  (GAZLC02),  8  ma3ffly  taxa  were  collected,  and  the  biotic 
index  value  (3. 17)  was  within  expectations  for  a  foothill  stream.  Four  cold-stenotherm 
taxa  were  present;  these  induded  the  sensitive  perlid  Doroneuria  sp.  and  the  caddisfly 
Dolophilodes  sp.  Cold,  clean  water  is  suggested  by  these  findings. 

Six  caddisfly  taxa  and  13  "dinger"  taxa  were  collected.  It  seems  unlikely  that 
sediment  deposition  appredably  limited  the  availability  of  stony  substrate  habitats. 
Overall  instream  habitat  quality  was  probably  good,  with  diverse  undisturbed  niches, 
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since  42  taxa  and  a  rich  predator  fauna  (9  taxa)  were  evident.  No  fewer  than  10  stonefly 
taxa  were  supported  at  the  site,  suggesting  that  reach-scale  habitat  features,  such  as 
streambanks,  riparian  zones,  and  channel  morphologic  elements  were  essentially 
intact.  Long-lived  taxa  were  well-represented  (5  tax^;  recent  catastrophic  dewatering  or 
sediment  pulses  can  probably  be  ruled  out.  All  expected  functional  elements  were 
present.  The  abundance  of  shredders  is  notable,  and  could  indicate  ample  riparian 
sources  of  large  organic  material  and  good  channel  retention  of  such  matter. 

Spring- fed  stream  systems  are  more  difBcult  to  evaluate  than  run-off  driven 
systems  since  neither  of  the  bioassessment  methods  in  ciurent  use  in  Montana  have 
been  tested  for  their  applicability  to  spring  streams.  Interpretation  of  individual  metrics 
is  difScult,  since  none  have  been  investigated  for  associations  with  habitat  features  or 
water  quality.  In  the  upper  reach  of  Blaine  Spring  Creek  (BLNSC03),  non-insects  and 
midges  dominated  the  benthic  assemblage;  most  abundant  among  these  was  the 
hydrobiid  snail  Amnicola  sp.  Because  of  unique  hydrologic  and  thermal  conditions, 
spring  streams  often  support  invertebrate  communities  with  low  proportions  of  insects, 
but  such  assemiblages  result  in  low  bioassessmient  scores.  The  animals  taken  at  this 
site  seem  to  represent  an  appropriate  favina  for  a  spring-fed  stream. 

Farther  downstream,  the  caddisfly  HeUcopsyche  boreaUs  dominated  the 
assemblage  at  site  BLNSCOl,  and  the  non-insect  fauna  was  limited  to  a  few  snails, 
clams,  and  worms.  Cool  water  temperatures  are  su^ested;  no  cold-stenotherms  were 
present.  "Clinger^  taxa  are  more  prevalent,  implying  larger  substrate  particle  sizes,  and 
clean  surfaces  not  appreciably  contaminated  by  fine  sediments. 

Spring-flow  influence  in  Blaine  Spring  Creek  apparently  persists  to  the  site  at 
Vamey  (BLNSC04),  where  the  amphipod  Gammarus  sp.  and  physid  snails  dominate  the 
benthic  assemblage.  The  absence  of  cold-stenotherms  and  common  occurrence  of  cool- 
to-warm  water  forms  suggest  a  thermal  regime  similar  to  the  upstream  sites.  Thirteen 
"dinger"  taxa  and  9  caddisfly  taxa  imply  clean  stony  substrates. 

CONCLUSIONS 

•  Water  temperatures  may  have  been  warmer  than  expected  at  the  sampled  site 
on  Buford  Creek  (BUFDCOl).  Nutrient  enrichment  was  possible  as  well. 

•  Both  sampled  sites  on  Elk  River  (ELKROl  and  ELKR02)  supported  assemblages 
suggesting  good  water  quality  and  diverse  instream  habitats. 

•  Antelope  Creek  (ANTLCOl)  yielded  a  diverse,  sensitive  benthic  assemblage 
characteristic  of  minimally  distxirbed  stream  sites. 

•  Excellent  water  quality  and  intact  instream  habitats  were  indicated  by  the 
invertebrates  collected  at  both  sites  on  Gazelle  Creek  (GAZLCOl  and  GAZLC02). 

•  Tested  methods  for  the  assessment  of  spring  streams  are  not  available.  When 
the  revised  index  (BoUman  1998a)  is  appUed  to  data  collected  at  the  3  sites  on 
Blaine  Spring  Creek  (BLNSC03,  BLNSCOl,  and  BLNSC04),  results  suggest 
moderate-to-severe  impairment.  However,  the  taxonomic  composition  of  the 
faunae  collected  at  these  sites  appears  to  be  typical  of  spring  streams  of  Western 
Montana.  Impairment  has  likely  been  overestimated  by  the  revised  metric 
battery.  The  Montana  DEQ  battery  (Biokantis  1998)  appears  to  have  given  the 
more  accurate  assessment  resvdts  for  Blaine  Spring  Creek  sites;  slight 
impairment  was  indicated  by  that  method. 

Habitat  Assessment* 

Figure  3  graphically  compares  total  habitat  assessment  scores  recorded  for  the  9 
sites  in  this  study.  Table  5  shows  the  habitat  parameters  evaluated,  parameter  scores 
and  overall  habitat  evaluations  for  the  sites  studied. 
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Fl^iue  3.  Total  habitat  assessment  scores  for  sites  on  the  Madison  River  drainage. 
August  -  September  2003 
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Overall  habitat  conditions  were  rated  optimal  or  sub-optimal  at  aU  sites  studied. 
At  the  Buford  Creek  site  (BUFDCOl),  sediment  deposition  was  noted  behind  boulders 
and  on  stream  edges.  However,  overall  habitat  condition  was  appraised  as  optimal. 

The  upper  site  on  Elk  River  (ELKROl)  received  optimal  or  sub-optimal  scores  for 
all  iastream,  streambank,  and  riparian  zone  parameters.  At  the  lower  site  (ELKR02), 
marginal  benthic  substrate  diversity  was  noted,  and  bank  stability,  bank  v^etation 
protection,  and  riparian  zone  width  were  judged  marginal  or  poor  on  the  left  side  of  the 
channel.  Overall  habitat  conditions  were  rated  sub-optimal  at  both  sites. 

Some  sediment  deposition  was  reported  at  cow^  crossings  and  obstructions  and 
constrictions  at  the  studied  site  on  Antelope  Creek  (ANTLCOl).  Marginal  streambank 
stability  was  assessed,  and  the  riparian  zone  was  noted  to  be  abbreviated  on  the  left 
side  of  the  stream.  Overall  conditions  were  rated  sub-optimal  here. 

Substrate  particles  were  emibedded  to  some  extent,  and  moderate  sediment 
deposition  was  observed  at  the  upper  site  on  Gazelle  Creek  (GAZLCOl).  Sub-optimal 
overall  conditions  were  reported.  At  the  lower  site  (GAZLC02),  aU  instream,  streambank, 
and  riparian  zone  parameters  were  rated  optimal  or  sub-optimal,  and  overall  conditions 
rated  optimal. 

Marginal  riffle  development  and  rather  monotonous  benthic  substrates  were 
noted  at  the  upper  site  on  Blaine  Spring  Creek  (BLNSC03).  Other  habitat  parameters 
received  high  scores,  however,  with  the  result  that  overall  habitat  conditions  rated  sub- 
optimal.  At  the  middle  site  (BLNSCOl),  aU  parameters  were  judged  optimal  or  sub- 
optimal,  and  overall  conditions  rated  optimal.  The  lower  site  (BLNSC04)  was  rated  sub- 
optimal.  There,  monotonous  benthic  substrates,  embeddedness,  and  sediment 
deposition  were  reported. 
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Habitat  assessment  vs.  bloassessment 

When  habitat  assessment  scores  are  plotted  against  bioassessment  scores,  the 
resulting  figure  provides  an  opportunity  to  evaluate  the  hypothetical  relationship 
between  habitat  int^rity  and  water  quality.  Both  factors  are  critical  and  interactive 
determinants  of  the  composition  and  functional  integrity  of  aquatic  invertebrate 
assemblages.  Presumably,  high  quality  habitat,  in  the  absence  of  impairments  to  water 
quality,  supports  functional,  diverse,  and  sensitive  invertebrate  assemblages;  these  are 
assemblages  that  attain  high  bioassessment  scores. 

Barbour  and  StriblLng  (1991)  have  h3TJothesized  that  diminishing  habitat  quality 
should  produce  predictable  diminishment  of  bioassessment  scores,  when  water  quality 
is  not  a  further  insult.  Figure  4  is  a  plot  of  habitat  assessment  scores  against 
bioassessment  scores  (revised  method)  for  the  sampled  assemblages  of  the  Madison 
River  drainage.  The  red  line  sup>erimposed  on  the  plot  represents  the  hypothetical 
relationship  between  habitat  quality  and  biotic  integrity  given  good  water  quality.  In  this 
model,  symbols  falling  in  the  upper  right  area  of  the  graph  would  represent  sites  veith 
high  scores  for  both  bioassessment  and  habitat  assessment;  according  to  this  model, 
these  wovdd  be  unimpaired  sites  both  in  terms  of  habitat  integrity  as  well  as  water 
quality. 

The  plot  in  Figure  4  places  5  studied  sites  in  the  Madison  River  drainage  in  this 
r^ion:  ELKROl,  ELKR02,  ANTLCOl,  GAZLCOl,  and  GAZLC02.  When  habitat  scores  are 
high,  but  bioassessment  scores  are  inordinately  low,  sites  fall  into  the  lower  right  hand 
area  of  the  plot.  According  to  the  model,  these  sites  support  invertebrate  assemblages 
that  are  impacted  mostly  by  impairment  to  water  quality.  The  plot  in  Figure  4  indicates 
that  4  Madison  River  drainage  sites  fall  into  this  area.  They  are  BUFDCOl,  BLNSC03, 
BLNSCOI,  and  BLNSC04.  Since  the  bioassessment  method  used  here  has  not  been 
demonstrated  to  be  accurate  for  spring  streams,  the  position  of  the  Blaine  Spring  Creek 
sites  (BLNSC03,  BLNSCOI,  and  BLNSC04)  cannot  be  reliably  interpreted. 

Plgoze  4.  Total  bioassessment  scores  plotted  against  habitat  assessment  scores  for 
sites  on  the  Madison  River  drainage.  August  -  September  2003.  (Barbour  and  Stribling 
1991). 
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